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Abstract 

By making use of some symmetry properties of the relevant Hamiltonian, 
two fundamental relations between the ferromagnetic magnetization and a 
spin correlation function are derived for the d{= 1, 2, 3)-dimensional Hub- 
bard model at finite temperatures. These can be viewed as a kind of Ward- 
Takahashi identities. The properties of the magnetization as a function of 
the applied field are discussed. The results thus obtained hold true for both 
repulsive and attractive on-site Coulomb interactions, and for arbitrary elec- 
tron fillings. 
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The Hubbard model was originally proposed for interpreting the origin of itinerant fer- 
romagnetism in transition metals within the framework of mean-field approxmations more 
than thirty years ago It is now widely thought to be the simplest model for correlated 

electron systems. This model is exactly solved using the Bethe ansatz in one dimension [|J, 
but the exact results for dimensions d > 1 are sparse (see, e.g. Ref. || for a review). 

Among those rigorous results favouring the existence of ferromagnetism, the Nagaoka- 
Thouless's celebrated theorem || is the first example to show the itinerant electron ferro- 
magnetism in a special limit (with an infinite repulsive Coulomb interaction and one hole) 
in the ground state. Later, Lieb |7J] presented another example for ferrimagnetism in this 
model with finite Coulomb repulsion at half-filling on asymmetrical bipartite lattices. Mielke 
and Tasaki || showed the existence of ferromagnetism in a class of Hubbard models with 
somewhat artificial flat-band systems on decorated lattices. Very recently Miiller-Hartmann 
H proposed an idea which may open a low density route towards the understanding of 
ferromagnetism in the Hubbard models on non-bipartite lattices. 

At finite temperatures, only few rigorous results are known so far. They seem to be 
disfavourable of ferromagnetism in this model. For instance, the absence of magnetic order- 
ings was rigorously proved in one and two dimensions [|1(|, while the possibility of magnetic 
long-range order (LRO) in the attractive Hubbard model was precluded in arbitrary dimen- 
sions [11], etc. As the rigorous results are rare in this case, more efforts should be devoted. 
On the other hand, since the Hubbard model possesses a few nontrivial symmetries, which 
could induce some physically important properties, it is quite interesting to explore them 
so that a better understanding on the underlying physics behind these symmetries in this 
model could be achieved. 

It is the purpose of this paper that, by means of some symmetry properties of the 
relevant Hamiltonian, two fundamental relations between the ferromagnetic magnetization 
and a spin correlation function are derived for the d(= 1, 2, 3)-dimensional Hubbard model 
at finite temperatures. These can be viewed as a kind of Ward-Takahashi identities. The 
properties of the magnetization as a function of the applied field are discussed. The results 
thus obtained hold true for both repulsive and attractive on-site Coulomb interactions, and 
for arbitrary electron fillings. A special solution is analyzed as an example. 

We start from the Hubbard model in the presence of an external applied field h (> 0) 
on a d(= 1, 2, 3)-dimensional lattice A with |A| sites. The Hamiltonian reads 

H=- J2 t xy (cl a c ya + H.c.) + J2 u xn x ^n xi -hJ2( n ^- n xi), (!) 

x,y(^A,a x£A xGA 

where c xa (c xa ) is the creation (annihilation) operator of an electron at site ieA with spin 
a(=t, I), U x is the on-site electron-electron interaction at x, and n xa = c X(T c xa is the number 
operator of electrons. The elements of the hopping matrix, t xy , satisfying t xy = t* x , survive 
for the short-range hoppings of electrons in practice. Note that t xx = 0. The total number 
of electrons of the system, say N = J2 x eAa n xa, is conserved. We have taken the Lande 
factor of electrons g=2, and = 1. The foregoing discussion is independent of the sign of 
U x , and sustains valid only for finite temperatures which implies that one can not extract 
any useful information for zero-temperature case from this study. 

Define the spin operators as S + = J2 x eA c l-\ c xi, S~ — (S + Y , and S z = \ J2x&A( n xi; — n xi)- 
They satisfy the usual SU(2) Lie algebra: S~] = 2S Z , and [S*^*] = tS ± . The 
ferromagnetic magnetization per site is defined as 
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m(h,T) = ±-(S*), (2) 

where (• ■ ■) = Tr[exp(— (3H) ■ ■ -]/Tr[exp(— /3H)] is the thermal average in the canonical 
ensemble, with the inverse temperature (3 = 1/T (ks = 1). 

Using the commutator [S^if] = ±2hS ± , and the property of cyclicity under the trace, 
we obtain the following basic equation: 

g(h } T)(l-e- 2(3h ) = 2m(h,T), (3) 

where Q(h,T) = (1/|A|) (S + S~)h,T is the spin correlation function. This holds generally 
when the Hamiltonian without a Zeeman term satisfies 0(3) symmetry in spin space. This 
relates an odd spin correlation such as the magnetization to an even spin correlation. Thus, 



this can be viewed as a kind of Ward-Takahashi identities |12| in quantum field theory. It 
should be noted that there are an infinite number of correlation identities [13[] in classical 
systems. It is evident that if we know the form of Q(h,T), then we can get m(h,T) from 
Eq.@. On the basis of this strategy, we now try to exploit such properties of the function 
Q(h, T) that it enables us to obtain an exact solution for Q(h, T). Invoking the spin up-down 
symmetry which exchanges the up spins and the down spins of electrons, 

c*T c *h 4t C U ( 4 ) 

we find H — > H = Hq + 2hS z and S + S~ — > S~S + , where we have denoted the Hamiltonian 
[Eq. ([!])] by Hq — 2hS z . An alternative closed form for this unitary transformation can be 



found in Ref . [0 . It is interesting to note that this transformation, unlike the usual particle- 
hole transformation, conserves the total number of electrons (while it conserves neither the 
number of electrons with up spins nor the number of electrons with down spins). It is this 
property that makes it possible for us to work in the canonical ensemble. By applying the 
transformation, Eq.([|), to Q(h,T), we obtain 

g(h,T) = G(-h,T)e 2 P\ (5) 

where we have utilized the cyclicity under the trace in Q(h,T), S~S + = S + S~ — 2S Z , and 
Eq.(^|) as well. This equation is of basic importance for discussing the poperties of the 
correlation function Q(h,T). 

Generally speaking, Eq. (^) has an infinite number of solutions, but the number of the 
physically meaningful solutions, which could depend on the details of the Hamiltonian of 
the system, might be limited. Considering this point, we shall in the following only search 
for those solutions which are physically sound. Obviously, it would not be possible to solve 
Eq. (||D without first obtaining the conditions satisfied by the function Q(h,T). For this 
purpose we differentiate Q(h,T) with respect to h, and using the unitary transformation (Hp 
to the thermal averages involved, and then noting Eqs. @ and ([5]), we obtain 

dQ(h, T) dg(-h,T) 

dh ~ dh ' 1 ' 

where we have utilized the facts [((^) 2 ) - (S z ) 2 } h , T = [({S z ) 2 ) - (S Z ) 2 }^ T and ({S z ) 2 ) - 
{S z ) 2 > 0. In addition, by Eqs. (||), (H) and @ we get the following inequalities: 
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*<s(h,T)< r -? Fm (7) 

and 

where p = N/\A\ is the density of electrons in the system. On the other hand, we can obtain, 
after performing the transformation (|j) to m(h,T), the equation 

m{h,T) = -m{-h,T), (9) 

which implies that m(h,T) must be an odd function of h at given T. 

Up to now we have had some basic information about the function Q(h,T), namely, 
Eqs.(|^) and (|5|) - (|9|) should be satisfied simultaneously. Under these conditions, we can 
solve Eq. ([5]) exactly. The result, being surprisingly simple, takes the following form: 

where the function f(h, p, T), whose properties are to be discussed below, is an even function 
of h. Substituting this solution into Eq. (|) one gets 

m(h, T) = /(/l, 2 P,T) tanh(/3/i). (11) 

Evidently, Eqs. ([10]) and (JTl|) satisfy Eq. (0) and Eq. (g), respectively. Eq.(0) should 
remain valid in the thermodynamic limit (N — > oo and |A| — > oo with the ratio N/\A\ = p 



dm(h,T) | 
dh 



fixed). As a consequence, the magnetic susceptibility can be obtained by x{T) 
Eq. (|TTT) forms the main result of this paper, which is somewhat general for the Hubbard 
model defined in Eq. ([[]). 

Now, we are in the position to discuss the properties of the function f(h,p,T). 

1. As mentioned before, f(h,p,T) is an even function of the applied field, namely, 

f(h,p,T) = f(-h,p,T). (12) 

2. From Eqs. (0) and ([|), one may find that for any values of h the function f(h,p,T) 
should satisfy the following two inequalities: 

< f{h,p,T) < pcoth(/3/i), (13) 

(Recall that h > 0.) 

3. From the physical point of view, we note that m — > m max = S to t/\M as h — > +oo 
at given T, because in this limit all spins of electrons should align on the direction of the 
applied field, where S to t is the total spin of the system, i.e., S to t = N/2 in our case. m max is 
therefore the saturation magnetization per site. On account of this consideration, we find 
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Urn f(h,p,T)=p. (15) 

n—*ao 

This constraint is a consequence of the well-known physical result: \m(h,T)\ < m max for 
any h and finite T. 

4. Siito [15j proved an asymptotic expression for the ferromagnetic magnetization per 
site in the Hubbard model: m(T, h, p) — ► tanh(/3/i) as p — ► 1 (note that there is a difference 
in the definition of the magnetization by a prefactor 1/2). By means of this result, we have 

hmf(h,p,T) = l. (16) 

5. The function f(h,p,T) might be dimension-dependent, i.e., it may probably have dif- 
ferent forms in different spatial dimensions. The reasons are as follows. Since the rigorous up- 
per bounds for m(h, T) have been known for small h by using Bogoliubov's inequality, namely 
\m(h,T)\ < (const/T 2 / 3 ) l/il 1 / 3 for one dimension, \m(h,T)\ < (const/T 1 / 2 ) (1/| In \h\ | 1/2 ) for 
two dimensions flIU|| , and \m(h,T)\ < (const/T 1 ^ 2 ) for three dimensions, and in addition, the 
magnetic susceptibility for U x < has also been bounded by a constant in arbitrary dimen- 
sions in Ref. JTTJ] , we observe that f(h — > + , p, T) could be finite or could be divergent but in 



latter case the speed of divergence must be slower than that of [tanh(/t — > + )] _1 — > oo in one 
and two dimensions, and in three dimensions for U x < as well, otherwise it contradicts the 
consequences of these upper bounds. On the other hand, as can be seen, there are still some 
rooms where f(h — > + ,p, T) might be divergent for U x > such that m(h — ► + ,T > 0) 
could sustain a finite nonzero constant in this situation. Although we have the condition 
(|T6j) which is valid for any dimensions and independent of U x , we still can not preclude such 
a possibility from the present formalism due to the p-dependence of f(h,p,T). 

Therefore, any solution for the function f(h,p,T) should comply simultaneously with 
the properties given in the items 1-5. Considering these conditions obeyed by the function 
f(h,p,T), one may check that the solutions given by Eqs. (jTOj) and flllD indeed satisfy all 
conditions imposed on the functions Q(h,T) and m(h,T) self-consistently. 

Now, let us present, as an example, a possible realization for the function f(h,p,T), 
which is a well-known result for limiting cases in the system. Such a solution for the function 
f(h,p,T) takes the simplest form, namely, 

f(h,p,T)=p, (17) 

which gives rise to the following ferromagnetic magnetization per site, m(h,T), in the Hub- 
bard model at finite temperatures: 

m{h,T) = ^tanh(/3/i). (18) 

It is easy to verify that this special solution satisfies all conditions imposed on the functions 
f(h,p,T), Q(h,T) and m(h,T) self-consistently. As can be easily seen, it is nothing but 
the result for the non-interacting system (U x = 0), as well as the result for the system at 



half-filling, as obtained by Siito |15| . It leads to the magnetic susceptibility obeying the 
Curie's law, and the fluctuations of the spin-spin correlation function in the presence of 
the applied field vanishing in the thermodynamic limit. We would like to point out here 
that we still can not determine if the form of m(h,T) given in Eq.([H|) pertains to the one- 
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and two-dimensional systems with arbitrary U x , and to the three-dimensional system with 
U x < from the present formalism, although we are aware of the fact that no spontaneous 
ferromagnetic LRO appears in these cases at finite temperatures. 
Several remarks are collected here. 

(i) As one may notice that, apart from the case at half-filling, we can not, only from the 



special solution given by Eq. (|lq) , draw the conclusion that no spontaneous ferromagnetic 



LRO emerges in the three-dimensional repulsive Hubbard model at finite temperatures, 
because there might be other possible solutions which could lead to the exhibition of the 
ferromagnetic LRO at finite temperatures in some parameter regimes. Since we use only the 
symmetry properties and known exact results but not invoking any detail information of the 
Hamiltonian, we can not at the moment affirm what solution is more physically realistic in 
the parameter regimes w hich are currently interesting. It is conceivable that resolving this 
problem the explicit form of the Hamiltonian ([!]) or numerical works ought to be involved 
in the analyses. However, we expect that our exact result [Eq. (|TT])] could shed some useful 
light on this issue. 

(ii) We should stress that since the Hubbard model can be mapped onto the t-J model 
in the limit of U x 3> t, while the t-J model at half-filling reduces to the antiferromagnetic 
Heisenberg model which may exhibit LRO at finite temperatures in dimensions d > 3, the 
result [e.g. Eq. (|ll]) or Eq.(|T8|) ] is not in conflict with this consensus. Actually, people can 
expect that the antiferromagnetic LRO would occur in the half-filled, repulsive Hubbard 
model at positive temperatures in three dimensions 0. 

(iii) We believe that the present result is suitable for one, two and three dimensions. It 
is still uncertain whether or not the method used in this paper can be extended to either the 



Hubbard model in infinite dimensions ||16|| , or the Hubbard models with artificially degener- 
ate bands on decorated lattices ||, or the spin systems which possess spin SU(2) symmetry, 
like the ferromagnetic Heisenberg model, and so on, because the additional symmetries or 
restrictions in these systems, if they exist, might violate this method in some ways. This 
question remains open. 

(iv) The present result holds even if we add an arbitrary one body potential (real), e.g. 
y = DxeA Vx( n xi + n xi), into the Hamiltonian [Eq. ([[])]. 

(v) We like to mention that although we can not prove the uniqueness of the solutions 
for Q(h, T) and m(h, T) at the moment, the other solutions, if they exist, might have similar 



forms as those given in this paper due to the restrictions imposed on the two functions []T7 



Despite all this, the present result at least is one of the exact results (if any) of the model. 
Moreover, our derivation, without invoking any a priori assumption, is somewhat general. 

(vi) It can be seen that any further solutions, exact or numerical, of m(h, T > 0) in the 
Hubbard model should comply with the fundamental equations and results discussed in this 
paper, otherwise they would conflict some symmetry properties of the model. 

In summary, by means of some symmetry properties of relevant Hamiltonian, two funda- 
mental relations between the ferromagnetic magnetization and a spin correlation function, 
which can be viewed as a kind of Ward-Takahashi identities, are derived, which allows us, 
combining some known exact results, to discuss symmetry properties for the ferromagnetic 
magnetization per site as a function of the applied field in the d(= 1, 2, 3)-dimensional Hub- 
bard model at finite temperatures. The results hold true for both repulsive and attractive 
on-site Coulomb interactions, and for arbitrary electron fillings. As an example a special 
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solution is also analyzed. The properties discussed in this paper could be helpful for further 
understanding the magnetic aspects of the Hubbard model at finite temperatures. 
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